%a=load('t1.txt');
clear all;
a=load('t3.txt');
c = ['r','b','k','g'];
%chr=['r','r','b','b','c','c','y','y','m','m'];

[row col]=size(a);
%t=[1:4];
%AXIS OFF;
%AXIS([-6000 12000 0 0.004]);
%AXIS([XMIN XMAX YMIN YMAX])
hold;
for i=1:2:col
%for i=3:2:3
%subplot(4,1,i)
subplot(col/2,1,(i+1)/2)
    for j=1:row
        mean=a(j,i); sigma=a(j,(i+1));
        Y=[(mean-3*sigma):0.1:(mean+3*sigma)];
        [P Q]=size(Y);
        mean=mean*ones(Q,1); mean=mean';
        sigma=sigma*ones(Q,1); sigma=sigma';
        X=[(1./(sigma*sqrt(2*pi)).*exp(-((Y-mean).^2./(sigma.^2))))];
        %plot(Y(1,:),X(1,:),'y-');
        %AXIS OFF;
        %AXIS([-2500 10500 0 0.0035]);
        plot(Y(1,:),X(1,:),c(1,j));
        %grid on
        hold on;
        %axis off;
        %r=a(j:j,:);
        %plot(r',a(:,i),'k');
        %subplot(4,1,)
        %clear X;clear Y;clear P;clear Q;
            
    end
   end
hold;




